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We propose phonon spectroscopy by electric measurements of the low-temperature conductance of coupled- 
quantum dots, specifically employing dephasing of the quantum electronic transport by the phonons. The setup 
we consider consists of a T-shaped double-quantum-dot (DQD) system in which only one of the dots (dot 
1) is connected to external leads and the other (dot 2) is coupled solely to the first one. For noninteracting 
electrons, the differential conductance of such a system vanishes at a voltage located in-between the energies of 
the bonding and the anti-bonding states, due to destructive interference. When electron-phonon (e-ph) on the 
DQD is invoked, we find that, at low temperatures, phonon emission taking place on dot 1 does not affect the 
interference, while phonon emission from dot 2 suppresses it. The amount of this suppression, as a function of 
the bias voltage, follows the effective e-ph coupling reflecting the phonon density of states and can be used for 
phonon spectroscopy. 

PACS numbers: 71.38.-k, 73.21.La, 73.23.-b 



I. INTRODUCTION 

Detecting dephasing sources in semiconductor quantum- 
dot devices, or alternatively investigating the hallmarks of 
quantum coherence, is of much importance for their various 
applications. The coherence of electrons passing through a 
quantum dot was demonstrated in a series of experiments, 1-3 
in which the dot has been embedded on an Aharonov- 
Bohm 4 (AB) interferometer. The wave of an electron travers- 
ing the arm of the interferometer carrying the quantum dot 
interferes with the wave passing through the other arm, re- 
sulting in AB oscillations in the conductance as a function 
of the magnetic flux penetrating the ring. Dephasing of the 
interference pattern in AB interferometers has been studied 
theoretically in several papers, in conjunction with electronic 
correlations, 5,6 due to coupling with an environmental bath, 7 ' 8 
or with phonons. 910 Another manifestation of coherence in 
AB interferometers was demonstrated in Ref. 3. It was ob- 
served that when high coherence is kept over the whole inter- 
ferometer, its conductance shows an asymmetric shape which 
is ascribed to a Fano resonance. 1 1 This resonance results from 
the interference of tunneling paths through the continuum of 
energy levels in the ring and the leads, with paths passing 
through the discrete states in the quantum dot. Upon in- 
creasing the bias voltage 3 the Fano resonance gradually takes 
the symmetric Lorentzian shape characterizing a Breit-Wigner 
resonance of the dot alone. Interactions of the transport elec- 
trons with phonons (while residing on the dot) have been 
shown to explain qualitatively the shape change of the Fano 
resonance with increasing bias voltages. 9 

Transport electrons passing through a quantum-dot system 
can emit and absorb phonons there. These electron-phonon (e- 
ph) interactions can be accompanied by energy exchange (in- 
elastic transitions) or not (i.e., when the same vibration modes 
are emitted and re-absorbed). These two types of processes 
play different roles in the transport properties of quantum-dot 
systems. The latter processes cause the "dressing" of the elec- 
trons, resulting in shifts of the energy levels in the dot. For 
example, these elastic processes narrow the resonance peak of 



the linear-response conductance plotted as a function of the 
gate voltage (this has been ascribed in Ref. LL2J to the Franck- 
Condon blockade, see also Ref. 13), and reduce the height of 
resonance peak at a finite bias voltage. In the case of inelastic 
processes, in particular at a finite bias voltage, the transport 
electrons emit or absorb phonons and consequently change 
their energy states. The inelastic processes diminish the life- 
time of the electrons on the dot, and hence contribute to the 
broadening of the conductance peak. At zero temperature, and 
for coupling with optical phonons, this requires a finite bias 
voltage. 

It is well-known that inelastic processes may lead to 
dephasing. 14 Consider for example an AB interferometer car- 
rying a quantum dot on one of its arms. When the electron 
emits (or absorbs) a phonon while residing on the dot, the in- 
terference between the waves passing by the dot and those 
which do not vanishes, because the corresponding phonon 
states are orthogonal to one another. However, inelastic pro- 
cesses do not always act as a dephasing source. When the in- 
terfering electron waves interact with the same inelastic scat- 
terer (emit or absorb the same phonon), that scatterer cannot 
be the origin of dephasing. This can be explained 14 using the 
example of the h/e and h/2e AB oscillations in the conduc- 
tance of AB rings. In the h/e case, one wave goes through one 
arm of the ring while the other wave goes through the other. 
When the scatterer is located on a certain arm, only one wave 
interacts with it. In the h/2e case, one wave goes around the 
entire ring and the other wave goes along the reverse direc- 
tion. Therefore, both waves interact with the same scatterer. 
In this case the states of the scatterer which are coupled to the 
waves are the same, and consequently the scatterer does not 
harm the interference. 

In this paper we examine a different type of a quantum- 
dot interfering device and show that there again one may en- 
counter a situation in which inelastic processes do not nec- 
essarily destroy coherence. When they do, though, one may 
exploit their effect to extract information on the effective cou- 
pling with the phonons, by measuring the differential conduc- 
tance as a function of the bias voltage. 
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FIG. 1: The T-shaped double-dot system. The upper dot, dot 1, is 
connected to the external leads by the tunneling matrix elements 
and t,R ,; dot 2 is connected, by the tunneling matrix element tc, solely 
to dot 1. Both dots are represented by a single energy level each, ei 
and £2. 



The setup we consider is the T-shaped double-quantum-dot 
(DQD) structure depicted schematically in Fig. [T] Represent- 
ing each quantum dot by a single energy level, e\ and £ 2 , one 
finds that when electron-phonon interactions are ignored the 
zero-temperature differential conductance of the device is 

dl 2e 2 aEi/4 

= =— = TT^- = • (1) 

dV h {E l E 2 -\t c \ 2 Y + 

Here, Ei = (e,; — eV)/T [i — 1,2] is the renormalized on- 
site energy on each of the dots relative to the bias voltage eV, 
measured in units of the level broadening due to the coupling 
with the leads, T = IY + Tr. The partial width resulting 
from the coupling with the left (right) lead is denoted by Tl 
(Tr). The asymmetry in these couplings is described by a — 
^lTr/F 2 , and tc = ic/I\ where tc is the tunneling matrix 
element coupling the two dots (see Fig. [TJ. Below, we choose 
the chemical potential on the left lead to be = eV, and that 
on the right one as /ir , = 0. 

The dotted line in Fig. [2] is the differential conductance 
computed from Eq. (Q]l for a symmetric DQD, e = e\ — 
£2 and a = 1. It has a double-peak structure with a dip in- 
between. The two peaks reflect resonant tunneling through the 
bonding and anti-bonding orbital states of the quantum dots, 
yielding dl/dV — 2e 2 /h. At the midpoint between the peaks 
the differential conductance vanishes due to fully destructive 
interference. This structure resembles a Fano resonance in a 
special casejii indeed the conductance may be fitted to the 
Fano line shape, (q + E) 2 j(E 2 + 7 2 ), where in our case, 
E = E 2 ,E 1 = 0,q = 0, and 7 = 2\t c \ 2 . 

In the presence of e-ph interactions, this double-peak struc- 
ture is modified significantly, as exemplified by the full curve 
in Fig. |2] Of particular interest is the effect of those inter- 
actions on the destructive interference leading to the dip: we 
find that only the e-ph interaction on dot 2 is responsible for 
the ascent of the dip. Moreover, we argue that the dependence 
of this ascent on the bias voltage reflects the effective cou- 
pling of the e-ph interaction on dot 2, and hence may serve for 
phonon spectroscopy. The effect of the e-ph interaction taking 



place on dot 1, together with the one on dot 2, is to decrease 
the peaks' height as compared to the zero-interaction case. 

Below, we study separately the case of e-ph interactions 
with acoustic phonons, and with optical ones. In semiconduc- 
tor quantum dots, the electrons interact with the bulk phonons 
of the semiconductor base. Then, the electronic interaction 
with acoustic phonons plays a major role at low tempera- 
tures since the energies of the optical phonons are much larger 
(~ 36meV for GaAs) than those of the acoustic phononsj^^ 
In the case of e-ph interaction with acoustic phonons the den- 
sity of available states is continuous, and indeed we find that 
the dip in the conductance rises up gradually as the bias volt- 
age increases. 

The interaction of transport electrons with optical modes is 
considered to be dominant in molecular junctions 17-21 which 
do not lie on a substrate. Several theoretical works study- 
ing e-ph interactions in molecular junctions have employed 
as a model of the molecular bridge a single-level quantum 
dot coupled to optical vibrations, calculating the differential 
conductanc e) 12 ' 13 ' 22 " — and the shot noise. 29-31 We find that the 
interaction with optical phonons affects the dip once the bias 
voltage matches the vibrational energy on dot 2, but not on 
dot 1. A delicate issue is the question of the population of 
the vibrational states. 28 When the e-ph interaction is with the 
acoustic phonons of the substrate, one may assume that those 
are given by the thermal-equilibrium distribution. This as- 
sumption does not necessarily hold for interactions on molec- 
ular bridges, where the population of the vibrational modes 
may be determined mainly by the transport electrons. Here 
we will assume that the molecular vibrations are thermalized 
by a coupling with the phonon bath of the surrounding. 

The organization of the paper is as follows. We begin in 
Sec. II by describing our model and presenting the expression 
for the differential conductance through the T-shape system. 
The technical details of the calculation are described in the 
Appendix, in particular the treatment of the electron-phonon 
interaction in the self-consistent Born approximation. 9 — For 
clarity, we confine ourselves to the case of zero temperature. 
Then, the transport electrons can only emit phonons once the 
bias voltage exceeds the phonon energy. Section III is devoted 
to the analysis of our results. First, we discuss the dephasing- 
free phonon emission. We show that the conductance dip is 
not affected by phonon emission from dot 1 . Second, we ex- 
amine the conductance dip as a function of the bias voltage 
V. We show that the amount of decrease in the dip as the bias 
voltage is increased follows the product of the phonons' den- 
sity of states and the coupling strength of e-ph interaction. We 
conclude by a summary/discussion section. 

II. THE MODEL AND THE CALCULATION METHOD 
A. The model 

Our model system is depicted in Fig. Q] and is described by 
the Hamiltonian 

H = H e \ + -He-ph + -Hph • (2) 
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with the acoustic phonon Hamiltonian being 



FIG. 2: The differential conductance dl/dV of a symmetric (e = 
£1 = £2, a — 1) DQD structure as a function of the gate voltage 
e. In the absence of the e-ph interaction, the conductance follows 
the dotted curve. The solid line is the differential conductance when 
the transport electrons are coupled to acoustic phonons on both dots. 
Here t c = 0.5T. 



Omitting the spin indices, the Hamiltonian of the electrons is 
H cl = H L + H R + H T + H D . (3) 
It consists of the leads' Hamiltonian 

H H&) = J2 £ fc C [(R)fc C L(R)fc: W 
k 

the DQD Hamiltonian 

H D = e 1 d\d 1 + e 2 d\d 2 , (5) 
and the tunneling Hamiltonian 

k k 

+ t c (d\d 2 +H.c). (6) 



Here, c 



and c T 



denote the creation and annihilation 



-L(R)fe ^L(R)fe 

operators of an electron of momentum k in lead L (R), respec- 
tively. We assume a single energy level E\ in dot 1 (e 2 in dot 
2), with the creation and annihilation operators on that level 
being d\, d 1 (d\, d 2 ). Electronic correlations are not included 
in our analysis. The broadening of the resonant level on dot 
1, T = Tl + Tr, is given by Fl(r) = where v is 

the density of states of the electrons in the leads, and (£r) 
is the tunneling matrix element connecting the dots to the left 
(right) lead. 

We consider e-ph interactions which are confined to the 
DQD. In the case of acoustic phonons, the interaction Hamil- 
tonian reads 

ffc-ph - M ?,iK + a - q ) d \ d i 

q 

+ ^M 9;2 ( ft? + fl g4d 2 , (7) 



(8) 



Here, a q and a q are the creation and annihilation operators 
of phonons with momentum q. We disregard the possibility 
of nonlocal e-ph interactions, (a q + a i _ q )dld 1 , assuming that 
the overlap between the wave functions of dot 1, and 
dot 2, \d 2 ), is small. The dispersion relation of the acoustic 
phonons is taken to be linear, 



cs\Q\ 



(9) 



with the sound velocity c$- When the e-ph interaction with 
acoustic phonons originates from the piezoelectric coupling, 
the matrix element of the coupling with the electron i s 15 ' 16 



M q>i = \ q {di\e iq - r \di 



with 
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(10) 



(11) 



For example, in GaAs g = 0.1. Because of the oscillating fac- 
tor e lq v , the e-ph interaction decreases once the wavelength 
of the phonons is smaller than the size of the dot Li(i = 1,2). 
For this reason one may choose 



\M, 



V2 



|A g | 2 
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TrVar? |q|2 



(12) 



(The effect of the mixed product M q i M* j for i ^ j is neg- 
ligible, because the phase-differnece between the two matrix 
elements that depends on the inter-dot distance diminishes its 
contribution.) 

An ubiquitous model to describe the interaction with optical 
phonons is the Frohlich Hamiltonian, 32 treating the quantum 
dots as Einstein oscillators of frequency Wj, i = 1,2. In that 
case the electron-phonon Hamiltonian reads 

H c - ph = Ci(a x + 4)4^1 + (2(02 + 4)4 d 2> (I 3 ) 
and the phonon Hamiltonian is 

-Hph = oo 1 a[a 1 + ui 2 a\a 2 , (14) 

where and a\ are the creation and the annihilation opera- 
tors of phonons on dot i and Q is the coupling strength of the 
optical e-ph interaction. 

Using the Keldysh formalism,— _ — the current can be ex- 
pressed in terms of the Fourier transform of the retarded Green 
function of dot 1, G r n (t-t') = -i9(t - t')U d i(*)> 4(*')}>> 

2 ' f du(~r)[f h (u)-f K (cj)}ImG r n (uj), (15) 



/ = 



see the Appendix for details. Here, /l(r)(w) = (exp[/3(w — 
A'l(r))] + 1) 1 is tne Fermi distribution function in lead L (R). 
In this paper we study the differential conductance, given by 



dG r lx (w) - 
' d(eV) 



(16) 
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where the explicit form of G^^ui), in terms of the self- energies, is (see the Appendix) 



£2 



(w - ei )(w - e 2 ) + |r(w - e a ) - \t c \ 2 - (w - E2 )£^(u;) - ^EL(^) ' 



(17) 



r 



m. RESULTS 

We begin with results pertaining to the case in which the 
electrons are coupled to acoustic phonons, see Eq. (0. In 
Fig. [3] we plot the variation of the differential conductance of 
a completely symmetric DQD with the gate voltage, scaling 
all energies by T, the level broadening on dot 1 . We present 
separately results for the case where the e-ph interaction takes 
place on dot 1 alone [panel (a)] or on dot 2 alone [panel (b)]. 
Both panels show also the differential conductance obtained in 
the absence of the coupling with the acoustic phonons. (The 
parameters chosen are given in the caption of Fig. 0) As is 



seen from Fig. [3] while phonon emission from dot 1 does not 
affect the dip in the differential conductance, it has a detrimen- 
tal effect on dl/dV when it occurs on dot 2. The destructive 
interference leading to the dip in the differential conductance 
is severely harmed by the e-ph interaction on dot 2. 

To explain this observation, we treat the e-ph interaction to 
second-order in perturbation theory (note, however, that the 
plots in Fig. [5] were computed in the self-consistent Born ap- 
proximation). The dip in dl/dV, Eq. ( TToT ). occurs at e = eV. 
Then, the first term there vanishes [since at zero temperature 
ImG^ (eV) — 0] while the second term yields 



dl _ 2e 2 aT 
dV ~ ~h~ 



eV 



lm( I tie 



( W _ eV f + -r( W - eV) \t c \ 2 - ( W - eV)X r u (u,) - -L2L_££ 3 ( W )_ 



deV 11 deV 

Calculating the self-energies appearing in this expression in second-order perturbation theory we find 



Eff»M=i;iM...i a /^i< ) ( 



{ 



LJ-U')\* 



Ld' — LUg + l0+ 



LJ' + LJ q + i0+ 



(18) 



(19) 



and 



uj' - uj q + i0+ 



LJ' + LJ q + i0+ 



(20) 



A straightforward calculation shows that <i£^ 2 ' (u>) / d(eV) at e = eV vanishes [this follows from Eq. (I17H . Hence, when the 
e-ph interaction on dot 2 vanishes, so does the differential conductance at the mid-point. On the other hand, 



e^ 2) H 

deV 



v \M q g\ 2 r L / 



2tt \t c \ 2 \lu - uj q - eV + i0+ lu + uj q - eV + i0+ 
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Tl / cs n 2 



uj-eV 



2(27r)V2 \t c \ 2y L 2 > (u-eV) 2 + (c s /L 2 ) 2 ' 
I 



(21) 



yielding that when phonon emission takes place on dot 2 the dip in the differential conductance is modified. 
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FIG. 3: The differential conductance dl/dV (solid lines) as a func- 
tion of the gate voltage e of a symmetric DQD (e = ei = S2, 
a = 1), for e-ph coupling with acoustic phonons. The bias volt- 
age is eV = 2T and tc = 0.5T (T is the level broadening on dot 1). 
Panel (a): the e-ph interaction is on dot 1 alone; panel (b): the e-ph 
interaction is on dot 2 alone. The dotted curves are the differential 
conductance when the e-ph interactions are absent. 



Since the e-ph interaction on dot 1 does not harm the inter- 
ference around e = eV, one may say that it does not cause de- 
phasing. One may monitor this (zero-temperature) dephasing- 
free phonon emission by varying the size of the DQD. As is 
mentioned above [see Eq. dl2t l. the efficacy of the e-ph in- 
teraction with acoustic phonons depends on the dot size: it 
decreases as the dot size increases. The effect of the two dots' 
sizes on the differential conductance is investigated in Fig. [4] 
In both panels, the solid curves pertain to the case of equal- 
size dots. The dotted curve in panel (a) shows the modification 
in the differential conductance brought about by increasing 
the size of dot 2 (thus making the phonon emission there less 
efficient). It is clearly seen that the decrease in dl/dV towards 
the dip is severely disturbed. On the other hand, when the size 
of dot 1 is increased [panel (b)] the differential conductance is 
almost unchanged. 

Next we consider the differential conductance in the case 
where the electrons are coupled to optical phonons. The 
main difference between this interaction and the coupling with 
acoustic phonons discussed above, is that now the transport 
electrons can emit real phonons and change their energies only 
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FIG. 4: The differential conductance dl/dV as a function of the 
gate voltage e of a symmetric DQD (s = ei = £2, a — 1), for e-ph 
coupling with acoustic phonons. The bias voltage is eV = 2F and 
tc = 0.5F (F is the level broadening on dot 1). Panel (a): changing 
the size of dot 2. The solid line is the differential conductance when 
both dot sizes are equal, L\ = L2 — cs/ (2. OF) and the dotted line 
is for the case L x = c s /(2.0F) and L 2 = c s /(0.5r). Panel (b): 
changing the size of dot 1. The solid line is the differential conduc- 
tance when both dot sizes are equal, L\ — L2 = cs/(2.0F) and the 
dotted line is f or the case L\ — cs/(0.5F) and L2 = cs/(2.0r). 



when the bias voltage exceeds the phonon frequency (at zero 
temperature). This is portrayed in Fig. [5] Switching-on the 
coupling to optical phonons when the vibration frequency on 
dot 2 is larger than the bias voltage almost makes no differ- 
ence in the shape of the differential conductance [see panel 
(a) of Fig. 0. The modifications in the peaks' heights are due 
to elastic processes taking place on both dots. On the other 
hand, when the bias voltage is large enough (as compared to 
the vibration frequency on dot 2), there is a significant effect 
on the dip in dl/dV, see panel (b) of Fig. [5] 

Finally we study the possibility to use the dephasing effect 
of the e-ph interactions on dot 2 for phonon spectroscopy. Fig- 
ure|6]depicts the dependence of the conductance at the dip on 
the bias voltage for coupling with acoustic phonons [panel (a)] 
and optical phonons [panel (b)]. The coupling of the electrons 
to acoustic phonons is characterized by the spectral function 
obtained from the product of their density of states and the 
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FIG. 5: The differential conductance dl/dV as a function of the 
gate voltage e of a symmetric DQD (e = e\ = £2, ct — 1), for 
e-ph coupling with optical phonons. Here tc = 0.5r (T is the level 
broadening on dot 1), and the strengths of e-ph coupling on the dots 
are Ci — C2 = 0.3F. S introduced in Eqs. dA14t and dA15t is 
S = 0.5r. Panel (a): the phonon energy in dot 1 (dot 2) is u)i = V 
(0J2 = 3r), and the bias voltage is eV = F. The dotted line indicates 
the conductance in the absence of e-ph interaction. Panel (b): the 
same as in (a) with eV = 3F. 



(acoustic) e-ph matrix element squared, 

xM = ^|M q , 2 | 2 ^-^) 



cs N 2 



(22) 
(23) 



Only the spectral function on dot 2 is presented, since (as was 
elaborated upon above) the conductance dip is affected by the 
e-ph interaction on that dot. The function x(w) is shown by 
the dotted curve of Fig. [6] Although it does not coincide with 
the conductance curve it does follow it, in particular at higher 
values of the voltage. 

Panel (b) of Fig. [6] shows the value of the conductance dip 
(as a function of the bias voltage) in the case where the elec- 
trons are coupled to optical phonons. The prominent feature 
here is the peak obtained when the bias voltage matches the 
vibration frequency on dot 2. The decrease in the conductance 
dip becomes more pronounced as the tunnel coupling between 
the two dots, tc, is decreased. This is because when this cou- 
pling is strong, separate effects of the dots on the conductance 
is more blurred. On the other hand, a smaller value of this cou- 




ev/r 



FIG. 6: The differential conductance dl/dV as a function of the 
bias voltage V for a symmetric DQD (e = ei = £2, a — 1), at 
the conductance dip (e — eV = 0). The tunnel coupling between 
the dots is tc = 0.3F (solid lines) and tc = 0.5T (dashed lines). 
Panel (a): coupling with acoustic phonons; the dotted line is x( w ) 
[Eq. d23H as a function of u (marked on the right-side vertical axis). 
Panel (b): coupling with optical phonons. The vibration energy in 
dot 1 is ui\ = r and on dot 2 is ui2 = 3F. 



pling enables the manifestation of the e-ph interaction on dot 
2 to become more distinguished. We do not plot the spectral 
function corresponding to the case of optical phonons since 
our calculation does not take into consideration the origin of 
the life time of these phonons. The broadening S introduced 
in Eqs. ( IA14t and (I A 1 5t is a free parameter, that in fact can be 
extracted by fitting the curves in Fig. |6](b) to the experimental 
data. 



IV. DISCUSSION 

In summary, we have studied the effect of electron-phonon 
interactions in a T-shaped DQD, employing the self-consistent 
Born approximation. The differential conductance of this de- 
vice is characterized by the appearance of a double-peak struc- 
ture due to interference. The dip in-between the two peaks is 
quite sensitive to the dephasing effect of phonon emission. We 
find that while phonon emission from dot 1 does not affect the 
conductance dip, phonon emission from dot 2 leads to its as- 
cent. Therefore, the electron-phonon interaction on dot 2, in 
particular the effective e-ph coupling given by the product of 
the phonons' density of states and the e-ph interaction, can be 
probed electrically, by monitoring the variation of the conduc- 
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tance dip. 

The conductance dip always appears at £2 = eV, even 
when the DQD is asymmetric (for example, e\ 7^ £2). Since 
the behavior of the self energy, which is responsible for the 
appearance of the dip, remains the same also in this case, our 
conclusions should hold for the asymmetric double-quantum- 
dot system. On the other hand, we expect that a finite tem- 
perature will smear the conductance dip (even in the absence 
of the electron-phonon interactions). When the temperature 
differs from zero, the behavior of the self energy is slightly (at 
low temperature) modified, and the distinction between the 
effects of the e-ph interactions on dot 1 and those on dot 2 is 
less clear. A finite temperature will also facilitate phonon ab- 
sorption processes. However, these alone are not expected to 
modify our main results considerably. One may hope though, 
that phonon spectroscopy via the monitoring of the conduc- 
tance dip is still possible at low enough temperatures. This 
in particular is so since we propose to monitor the differential 
conductance at finite bias voltages. 

An interesting point is the experimental possibility to have 
the electron-phonon interaction taking place mainly on one 
of the dots forming the DQD system. A phonon bath may 
be realized by an electronically insulating hard substrate. It 
is harder to imagine how a quantum dot can be effectively 
detached from a phonon bath and become "floating"; this 
seems to be more realistic in the case of a molecular junction, 
for which, under these circumstances, the relevant vibrational 
modes are optical. Then, a sufficiently high phonon frequency 
(higher than the bias voltage required for the dip) will satisfy 
the above requirement. 
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Appendix A: Details of the calculation 

Our calculation requires all four Keldysh Green 
functions^— belonging to the DQD: the time-ordered 
one 

G t ij {t-t') = -i{Td i {t)d]{t')), (Al) 
the anti time-ordered one, 

G\At-t') = -i{td i {t)d\{t')), (A2) 



with T (T) being the time-ordering (anti time-ordering) oper- 
ator, and the lesser and greater Green functions, 



G< j (t-t') = i{d](t')d i (t)), 
G>(t-t') = - l {d l (t)d]{t')). 



(A3) 



Here i and j take the values 1 and 2. 

It is convenient to present the Fourier transforms of these 
Green functions in a matrix form 



(A4) 



where 7 = t, t, <, or >. The retarded and advanced Green 
functions follow from these functions, 



G r {uj) = G\lu) - G<(lo). 



(A5) 



All the necessary Green functions are obtained from the cor- 
responding Dyson equations, in which there appears the self- 
energy due to the e-ph interaction, E, 

G r (u;) = G r(0) + G r(0) (a;)S r (w)G r (w), (A6) 

and the lesser Green function of the DQD 

G<(w) = G r (w)£ < (w)G Q (w) 

+ [1 + G r (u;)£ r H]G <(0) (u;)[l + S Q (o;)G a (w)], (A7) 

where G^ is the Green function in the absence of the cou- 
pling to the phonons. 

We treat the e-ph interactions in the self-consistent Born 
approximation. 910 Since we focus on the effect of the inelas- 
tic electron-phonon processes, we discard the Hartree term. 
When the transport electrons are coupled to acoustic phonons, 
the required self-energies are 

S « = 7-Z>/,,i| 2 / doo'GUuj-uj^D'iq^'), (A8) 



and 



-]>>^| 2 / du'G<{u-u')D<{q,u>). (A9) 



Here D denotes the Fourier transform of the zero-order (i.e., 
in the absence of the coupling to the electrons) phonon Green 
functions, 



D\q,uj) = -2m[N q S(u> + u> q ) + N q 5(u) - oj q )] 
1 1 

UJ — LU q + i0 + U) + LU q — i0 + ' 



(A10) 



and 



D<(q,u) = -27Ti[{N q + l)6{oj + oj q ) 

+ N q S{u> - ui q )], (All) 

where N q is the phonon population of mode q. It is implic- 
itly assumed that the vibrational modes are equilibrated by 
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the coupling to another heat bath, such that their population is 
given by Bose-Einstein distribution function, l/[exp(/3w q ) — 
1]. Then, N q = at zero temperature. The self-energies 
Yj r u {oj) and E*^ are determined by solving Eqs. (IA8b and (|A9t 
self-consistently. 

When the electrons are coupled to optical phonons, the self- 
energies are 



yt _ J_j-2 

11 2tt 2 



du'GUuj-^DKq^'), 



and 



2tt 



C 2 / du'G<{u-u>')Df(q,u'l 



(A12) 



(A13) 



where the Fourier transforms of phonon Green functions at 
zero temperature are 



i 



i 



UJ — UJi + i S 



i5' 



and 



+ cjj - id 



(A14) 



(A15) 



Here S is relaxation rate of the Einstein phonon mode due to 
the coupling with the surrounding bulk phonons. As in the 
acoustic-phonon case, S^(a;) and EJj are determined by solv- 
ing Eqs. (IA12I ) and dAl 3b self-consistently. 

The Green functions of the DQD, in particular G\ x , deter- 
mine the expression for the current. The operator of the cur- 
rent between the left lead and the quantum dots is given by the 
time derivative of the electron number operator in that lead, 
N ^ = E fc c[ k c Lk , 

h = -2e(N L ) = -2ie([H,N L }) = 4eRe^t L G< Lfc (i,i) 

k 

4e, 



-Re 



dwi L ^G< Lfe (u;) 

k 



(A16) 



adding a factor of 2 for the spin components. Here Gf Lk (t — 

t') = i(cLfc(i')^i(*)} i s me le sser Green function and 
Gf Lk (ui) is its Fourier transform. The current from lead R 
to the quantum dots, J R , is obtained in the same way. 

Using the equation-of-motion method 35 - the current from 
lead L (R) is rewritten as 



4(R) =-£ I <!■*-• 



L(R) 



- r L (R)/L(w)ImGn(w) 



(A17) 



where G r denotes the retarded Green function. The net cur- 
rent through the DQD system is hence 



II - Ir 



T R ImG<H 



±(rMu) - r R / R H)lmG^H j. (A18) 



Using the relation Jl + Ir = which follows by charge con- 
servation, and assuming that the widths IY and T R do not 
vary significantly with energy, we can express the lesser Green 
function in terms of the retarded one, 



ImGf^Lj) = - 



2[r L / L n + r R / R H] 
r L + r R 



ImG^(w). (A19) 



Inserting Eq. JA19t in Eq. ( IA18b leads to Eq. ( fT3T l. 
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